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We investigate the vacuum realignment for chiral symmetry breaking and color superconductivity 
at finite density in Nambu-Jona-Lasinio model in a variational method. The treatment allows us to 
investigate simultaneous formation of condensates in quark antiquark as well as in diquark channels. 
The methodology involves an explicit construction of a variational ground state and minimisation of 
the thermodynamic potential. Color and electric charge neutrality conditions are imposed through 
introduction of appropriate chemical potentials. Color and flavor dependent condensate functions 
C^*) ' are determined through minimisation of the thermodynamic potential. The equation of state is 

calculated. Simultaneous existence of a mass gap and superconducting gap is seen in a small 
■ window of quark chemical potential within the model when charge neutrality conditions are not 

imposed. Enforcing color and electric charge neutrality conditions gives rise to existence of gapless 
| superconducting modes depending upon the magnitude of the gap and the difference of the chemical 

Q . potentials of the condensing quarks. 
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IT) ; 1. INTRODUCTION 

o i 

The structure of vacuum in Quantum Chromodynamics (QCD) is one of the most interesting questions in strong 
interaction physics [l| . The evidence for quark and gluon condensates in vacuum is a reflection of its complex nature 
■ 0> whereas chiral symmetry breaking is an essential feature in the description of the low mass hadron properties. 
Due to the nonperturbative nature of QCD in this regime different effective models have been used to understand 
the nature of chiral symmetry breaking ||. These have been constructed, for the most part, in the framework of a 
(-} ■ Nambu-Jona- Lasinio (NJL) model with a four fermion interaction. 

Recently there has been a lot of interest in strongly interacting matter at high densities. In particular, a color 
superconducting phase for it involving diquark condensates has been considered with a gap of about 100 MeV . The 
studies have been done with an effective four fermion interaction between quarks Q, direct instanton approach 
or a perturbative QCD calculation at finite density There has also been a study of this phase in NJL model 
Possibility of diquark condensates alongwith quark antiquark condensates has been considered in Ref. [M&EIEIE 
The natural place to look for such a phase seems to be in the interior of compact stellar objects like neutron star. 
However, to apply it to the case of neutron stars, the color and electric charge neutrality conditions need to be imposed 
for the bulk quark matter. Such an attempt has been made in Ref. 01 as well as i n Ref- [3 where the lighter up 
and down quarks form two flavor color superconducting (2SC) matter while the strange quark do not participate in 
pairing. It has been shown, based upon comparison of free energy that a two flavor color superconducting phase would 
be absent in the core of neutron stars 0] . Within NJL model in Ref. [l4| it has been argued that such conclusions are 
consistent except for a small window in density ra nge where superconducting phase is possible. There have also been 
studies to include the possibility of mixed phases 15] of superconducting matter demanding neutral matter on the 
average. In this context, there have been attempts to study the implications of vector interactions on the structure 
of the phase diagram [l^l ■ 

We had applied a different approach to study the problem in Ref. 9] . We considered a variational approach with 
an explicit assumption for the ground state having both quark antiquark and diquark condensates. The actual 
calculations are carried out for the NJL model such that the minimisation of the free energy density determines which 
condensate will exist at what density. In the present work we generalise the approach of Ref. 9] to include the 
conditions of color and electric charge neutrality. This leads to condensate functions which depend upon both color 
and flavor. Although, for simplicity, we shall be considering color superconductivity, this can be generalised to the 
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three flavor case to include color flavor locking. In fact here we shall also consider three flavors but having the u and 
d quarks taking part in diquark condensation. Although it might look rather complicated at the outset, the niceity 
of the approach is that within the model one can solve for the condensate functions explicitly which are flavor and 
color dependent. 

We organize the paper as follows. In the next section we discuss the ansatz state with quark antiquark as well as 
diquark condensates. In section 3 we consider Nambu Jona-Lasinio model Hamiltonian and calculate the expectation 
value with respect to the ansatz state to compute the thermodynamic potential. We minimise the thermodynamic 
potential to calculate all the ansatz functions and the resulting mass as well as superconducting gap equations here. 
In section 4 we give the results and discuss them. Finally we summarise and conclude in section 5. 



2. AN ANSATZ FOR THE GROUND STATE 



As noted earlier we shall include here the effects of both chiral symmetry breaking as well as diquark pairing. For 
the consideration of chiral symmetry breaking, we denote the perturbative vacuum state with chiral symmetry as |0). 
We shall then assume a specific vacuum realignment which breaks chiral symmetry because of interaction. 

Let us note first the quark field operator expansion in momentum space given as |0, ^} 



V>(x) 
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The superscript indicates that the operators q® and q® are two component ones which annihilate or create quanta 
acting upon the perturbative or the chiral vacuum |0). We have suppressed here the color and flavor indices of the 
quark field operators. The function <fi (k) in the spinors in Eq.(0l are given as cotcj)® = m;/|k|, for free massive 
fermion fields, i being the flavor index. For massless fields 0°(|k|) = ir/2. 

We now consider vacuum destabilisation leading to chiral symmetry breaking [tj 0, described by, 



U Q \0), 



where 



U Q = cxp (J qf (k)t( CT ■ k)hi{k)qf{k)dk - h.c. 
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(4) 



In the above, /i;(k) is a real function of |k| which describes vacuum realignment for quarks of a given flavor i. We 
shall take the condensate function h(k) to be the same for u and d quarks and hs^k) as the chiral condensate function 
for the s-quark. Clearly, a nontrivial /ij(k) shall break chiral symmetry. Sum over three colors and three flavors is 
understood in the exponent of Uq in Eq. (QJ . 

Having defined the state as in Eq.Q for chiral symmetry breaking, we shall next define the state involving diquarks. 
We note that as per BCS result such a state will be dynamically favored if there is an attractive interaction between 
the quarks |18|. Such an interaction exists in QCD in the qq color antitriplet, Lorentz scalar and isospin singlet 
channel. In the flavor antisymmetric channel the interaction can be scalar, pseudoscalar or vector whereas in flavor 
symmetric channel only the axial vector channel is attractive. In the present work, we shall consider the ansatz state 
involving diquarks as 



where 



|fi) = U d \vac) =cxp{Bl- B d )\i 



l (k)tr/-(k)^ r (-k)t £y - e3afc + ^(k)r/i°(k)e 6 r (-k)e ije3 a. 



dk. 
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In the above, i,j are flavor indices, a, b are the color indices and r(= ±1/2) is the spin index. As noted earlier we 
shall have u,d (i=l,2) quark condensation. We have also introduced here (color, flavor dependent) functions f ta (k) 
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and /i°(k) respectively for the diquark and diantiquark channels. As may be noted the state constructed in Eq.JSJ) is 
spin singlet and is antisymmetric in color and flavor. Clearly, by construction f la (k) = / J (k) with i ^ j and a ^ b. 
The corresponding Bogoliubov transformation for the operators is given by 



C'(k) 
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2re fcl e 3ca sin/ M (k) cos/ fcc (k) 



C(k) 
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In a similar manner one can write down the Bogoliubov transformation for the antiquark operators corresponding to 
|Q) basis. 

Finally, to include the effect of temperature and density we next write down the state at finite temperature and 
density |fi(/3, //)} taking a thermal Bogoliubov transformation over the state |fl) using thermofield dynamics (TFD) 
as described in ref.s [l|j|2(|- We then have, 



where Up^ is 



with, 
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In Eq. (|10|l the ansatz functions 9± (k, /?, /x) will be related to quark and antiquark distributions and the underlined 
operators are the operators in the extended Hilbert space associated with thermal doubling in TFD method. In 
Eci. (|10fl we have suppressed the color and flavor indices on the quarks as well as the functions 9(k,(3,fi). Note that 
we have a proliferation of functions in the ansatz state |0, /3, ji) - the (flavor dependent) chiral condensate function, 
the (color flavor dependent) quark as well as antiquark condensate functions and the (color flavor dependent) thermal 
functions. All these functions are to be obtained by minimising the thermodynamic potential. This will involve an 
assumption about the effective Hamiltonian. We shall carry out this minimisation in the next section. 



3. MINIMISATION OF THERMODYNAMIC POTENTIAL AND GAP EQUATIONS 

We shall work here in a Nambu-JonaLasinio model which is based on relativistic fermions interacting through local 
current- current couplings assuming that gluonic degrees of freedom can be frozen into point like effective interactions 
between the quarks. The Hamiltonian is given as 



n = ^ iot H" ■ v + i°mW a + Y J Z JfJ,a - ( n ) 

z, a 

Here J° = ipj^T a il) and, rrii is the current quark mass which we shall take to be nonzero only for the case of strange 
quarks («=3). Of the two superscripts on the quark operators the first index T refers to the flavor index and the 
second index, a refers to the color index. The point interaction produces short distance singularities and to regulate 
the integrals we shall restrict the phase space inside the sphere |p| < A the ultraviolet cutoff in the NJL model. 

We next write down the expectation values of various operators in the variational ansatz state given in Eq. JSJ. 
Using the fact that the state in Eq. © arises from successive Bogoliubov transformations one can calculate these 
expectation values. These expressions would be used to calculate thermal expectation value of the Hamiltonian to 
compute the thermodynamic potential. With -0(k)as defined in Eq.QJ, we evaluate the expectation values 

TO,M)IC°(k)<(k0 + |O(^ (12) 

and, 
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where, 



A£(k,/3,M) = 2 [l± (^r(k) -F iQ (k)) ± ( 7 °cos0 l (k)) + a • ksin^(k)) (l - F w (k) - ^(k)) 



(14) 



Here, the effect of diquark condensates and their temperature and/or density dependences are encoded in the functions 
F ia (k) and Ff°(k) given as 



and, 



F ta (k) = sin 2 6> w (k) + sin 2 f ia (Cl a (k) - sin 2 6>!_ a (k)) (l - 5 a3 ) 



F{ a (k) = sin 2 6^ (k) + sin 2 /*» (cos 2 0£(k) - S; a (k)) (l - 5 a ~ 



(15) 
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Here, we have defined C!_ a = |e M e aa | cos 2 l _ a and S^ a = \e vl ' t aa ' \sm z l + a ' . The 6 aa term indicates that the third 



e"'e oa '|cos 2 6» l >' and = \e vl 'e aa '\ sin 2 d% a ' . The c5 a3 
color does not take part in diquark condensation. Further, we have introduced the notation 0i(k) = </>°(k) — 2/i;(k). 
We also have 



(fi(/3,M)|C a (xX(0)|O(/3, M )) 
(^(/3, M )|C at (x)V4 fct (0)|^(/3,M)> 



1 



where, 



n ajb (k,/?, M ) 



ij ~3ab 



e' J e 



7° cos 



S ia ' jb (k)< 



(2tt) 3 
1 

(2^)3 



2 

a ■ k sin 



J e * x ^(k,/3,/i)dk, 



(17) 



A mjb (k) 



7 5 C, 



(18) 



and, 
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Here, C — i^ 2r y° is the charge conjugation matrix (we use the notation of Bjorken and Drell) and the functions S(k) 
and A(k) are given as, 



S lJ ' ab (k) = sin2/ w (k) cos 2<f a jb (k, 0, fi) + sin2/f(k) cos 29™ jb (k, /3, fi), 



and, 



A l3 - ab {k) = sin2/ M (k)cos20!_ ajb (k,/^) - sin2/j a (k) cos 26^ aj0 (k, 0, /x) 
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In the above we have defined cos 29±'-' b = 1 — sin 2 6± — sin 2 #£ b , with i, j = 1,2 being the flavor indices and a,b = 1,2 
being the color indices and i 7^ j, a 7^ b. 

Using Eq. I|13|) we have for the kinetic energy of the light quarks 

= T^Ta" E / ^k|k| (l — cos2/ij(k)(l — F ia — -F^ )) , (22) 

where, F la and F™ are defined in equations (|15|) and Ea. i|lt)|) . We have also subtracted the vacuum contributions. 

Similarly the contribution from the interaction term in Ea. (|ll|) after subtracting out the zero point perturbative 
energy turns out to be 



v = (n(0,fx)\^-j^ a \n(0, t ,)) = v 1 + V2 



(23) 



5 



Here, the contribution V\ arises from contracting -0 with a fusing Eq.s (|12|l . (|13() and is given as 



* = tE ( Ec 2 -2^/r 2 ), (24) 

i=l,2 \a=l,3 o=l,3 / 

with, 

C = (2^W / dk{Fta ~ ¥{a) (25) 

and 

IT = 7A3 | - - sin 2^(k). (26) 

In this expression, we have neglected terms of the order 1/N 2 compared to unity. The term V2 arises from contracting 
i/j and a ^ and i/^ with another ijft using Eq. Q18|l and (|19fl and we have 

^ = -| ff 2 /3 U ' 22 / 3 12 ' 21 (27) 

with 



r ia,jb 



(2tt) 3 



dk5 ia ^(k)cos( (28) 



where, 5" iajb has been defined in Ea. l|2U|l . 

To calculate the thermodynamic potential we shall have to specify the chemical potentials relevant for the system. 
Here we shall be interested in the form of quark matter that might be present in compact stars older than few 
minutes so that chemical equilibriation under weak interaction is there. The relevant chemical potentials in this case 
then are the baryon chemical potential n B = 3//, the chemical potential fiE associated with electromagnetic charge 
Q = diag(2/3, —1/3, —1/3) in flavor space, and the two color electrostatic chemical potentials /13 and (1$ corresponding 
to 1/(1)3 x U(l)s subgroup of the color gauge symmetry generated by cartan subalgebra Q3 = diag(l/2, —1/2, 0) and 
Qs = diag(l/3, 1/3,-2/3) in the color space. Thus the chemical potential is a diagonal matrix in color and flavor 
space, and is given by 

faj,ab = (M<% + QijHE)Sab + (Qzab + Qsab^^ij- (29) 

Here, i, j are flavor indices and a, b are color indices. 

The total thermodynamic potential, including the contribution from the electrons, is then given by 

n = T + V - (fiN) - -s + n e (30) 

P 

where, we have denoted 

(jxN) = (V wt A%-,a^ b ) = 2 £ » iai v ( 31 ) 

i,a 

with a la being the chemical potential for the quark of flavor % and color a, which can be expressed in terms of the 
chemical potentials fi, fiE,H3 and /xg using Eq.J55J. fl e = — u%/12tt 2 is the electron free energy. 
Finally, for the entropy density for the quarks we have |19| 



Yy Yl J dk ( sin2 e - ln sin2 e - + cos2 e - ln cos2 e - 



(2 

+ sin 2 6^ In sin 2 6^ + cos 2 0™ In cos 2 9™ ) . (32) 
Now if we minimise the thermodynamic potential f2 with respect to /ii(k), we get 



(33) 
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where, Mj = m; + g 2 J^ fc 7J b and e; = y (k 2 + mf). Substituting this back in Eg. 126(1 we have the mass gap equation 
for the light quarks (rrii—0) 



(34) 



Clearly, the above includes the effect of diquark condensates as well as temperature and density through the functions 
F and F\ given in Eq.s l|15f) and (|16|l respectively. 

Next, minimising the thermodynamic potential with respect to the diquark condensate functions leads, to 
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In the above e = (ej + ea)/2, = (^n + vii)/2, &12 = (^12 + ^2i)/2- is the interacting chemical potential given 
as 

„2 



v = (1 



-P 



(37) 



which may be expected in presence of a vectorial current-current interaction . In the above, p 1 = 2 X) a =i 3 W1 tfi Z* a 
as defined in equantion 125[) . Thus, it may be noted that the diquark condensate functions depend upon the average 
energy and the average chemical potential of the quarks that condense. We also note here that the diquark condensate 
functions depends upon the masses of the two quarks which condense through the function cos ((fa — fa)/ 2) with 
cosfa = sin2/ii = Mi/ti, for u,d quarks which could be different when charge neutrality condition is imposed. Such 
a normalisation factor is always there when the condensing fermions have different masses as has been noted in Ref. 
[2l| in the context of CFL phase. 
In an identical manner the di-antiquark condensate functions are calculated to be 
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Further, in Eq.s J33> , J3SJ , ISHl , J23) , A12 = (2# 2 /3)/ 3 u ' , An = (2 5 2 /3)J 3 12 ' 21 which satify the equations 
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In the above, £,±i a = e ± v %a . Finally, minimisation of the thermodynamic potential with respect to the thermal 
functions 0±(k) gives 
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exp(/3u/3?) + 1 



(42) 



7 



Various co la s are given as follows. 
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and, finally, for the noncondensing colors ui 

second index refers to color. Here u>± = J A 2 cos 2 ((^i — (f>2)/2) + 5 e = (e 



As already mentioned, the first index refers to flavor and the 
£2)/2 is half the energy difference of 

the two quarks which condense and e.g. 5* 1 = (v\\ — 1^22) /2, is half the difference of the chemical potentials of the 
two quarks which condense. Note that in the absence of imposing the charge neutrality condition all the four quasi 
particles will have the same energy cj_. Thus it is possible to have zero modes when charge neutrality condition is 
imposed depending upon the values of 5 e and S u . So, although wc shall have nonzero order parameter A, there will 
be fermionic zero modes or the gapless superconducting phase [23|,|^J. We shall discuss more about it section 4. 

Next, let us focus our attention for the specific case of superconducting phase and the chemical potential associated 
with it. First let us note that the diquark condensate functions depend upon the average of the chemical potentials 
of the quarks that condense. Since this is independent of /Z3 we can choose /L/,3 to be zero. In that case, fi 11 — 



2/3fiE + M8/V3 = fi 12 and also (j? 1 — fi — 1/3/is + [is/VS 



22 



This also means that the average chemical 



potential of both the condensing quarks are the same and is equal to ji = /1 + 1/6/iE + fAs/VS. Further 5 U 12 = 
\xeI^— (s ,2 /4)(/Oi — P2)/2 = 5 U = 8 V . In such a situation, Ea. (|4U|l and Eg. (|41|l are both identical and hence we shall 
have only one superconducting gap equation given as 
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With this condition, it is clear from Ea.s l|43|l that the quasi particle energies for each flavor becomes degenerate 
for both the colors which take part in condensation. Thus the quasi particle energies now become u\ = uj- + 8 and 
uj2 = W- — S, for u and d quarks respectively, with S = S e — <5„. 

Now making the use of this dispersion relations and the mass gap equations Ea. H34|) and the superconducting gap 
equation Eq. l44|) . the thermodynamic potential at zero temperature becomes, with 5 = 8 t — 5 V , 
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The first three lines in Eq. (|45|l correspond to the contribution from the quarks taking part in the condensation while 
the last line is the contribution from the third color for the two light quarks. 

In an identical manner one can calculate the thermodynamic potential for the strange quark sector with strange 
quarks and anti-quarks for the vacuum structure. The contribution to the thermodynamic potential then is given by 



(27T) 



a=l,3 
2 „2 



\k\ 2 + m 2 - Jk 2 + Aff 0(|k| - kf ) 



(M 3 - m s ) 2 g 
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Here, Pl = £ a=lj3 p M and P™ = 2I T for ^ a = with P v a given in Eq.|23). For the third color, p a = {kf) 3 /3n 2 
and for the strange quark, p 3a — (k 3a ) 3 /3n 2 . The fermi momenta are given by the usual relation k™ = [y 2 a — M 2 ) 1 / 2 . 
The total thermodynamic potential is given by including the contribution from the electrons and is given by 

n = n u , d + n s --^-. (47) 

Thus the thermodynamic potential is a function of four parameters: the three mass gaps and a superconducting gap 
which needs to be minimised subjected to the conditions of electrical and color neutrality. The electric and charge 
neutrality constraints are given respectively as 

2 11 

Qe = 3P1 - 3P2 - g/?3 - p e = 0, (48) 

and, 

Qs = ^J2(p tl + P l2 -2p a ) = 0. (49) 

Eo . (|47tf4*9")l and the superconducting gap equation Eq. l|44() constitute the basis of the numerical calculations that we 
discuss below. 



4. RESULTS AND DISCUSSIONS 



For numerical calculations we have taken the values of the parameters of NJL model as follows: g 2 A 2 = 17.6, 
A = 0.68 GeV as typical values giving reasonable vacuum properties 0, ■ We might note here that the 

coupling g 2 introduced here is related to the usual scalar coupling of NJL model e.g. in Ref.|22] G as g 2 — 8G. 

Current quark masses for u and d quarks are taken as zero and the current quark mass for strange quark is taken 
as 0.12 GeV. With this choice of parameters, the constituent quark masses at zero temperature and density are given 
as Mi = 0.35 GeV=M 2 , and for strange quark M 3 = 0.575 GeV. 

In Fig^ we have plotted the variation of masses with baryon density without diquark condensation and without 
imposition of the charge neutrality conditions. The decrease of masses with density is a reflection of the decrease of 
quark condensates with density. 

Let us next discuss the case with diquark condensates along with quark antiquark condensates without imposing 
the charge neutrality condition. The numerical calculation for this case proceeds as follows. For a given chemical 
baryon potential, (and with /i3=0= ps), the thermodynamic potential Ea. l|47|) is minimised with respect to the quark 
masses subjected to self consistently determining the interacting chemical potential using Eq. l|37|l and solving the 
superconducting gap equation Eg . (|44|l . In Fig. 2 we have plotted the masses and the superconducting gap as 
a function of baryon chemical potential. The behaviour of masses are almost same to that without the diquark 
condensates. Here, however, we observe that for a small window in the quark chemical potential ( about 23 MeV), 
we have simultaneous existence of chiral symmetry breaking and color superconducting phase with the gap reaching 
a maximum upto 65 MeV in this window. Including vector interactions similar conclusion of simultaneous existence 
of both the condenstes was noted earlier in Ref . [l2T| . However, in Ref.01 the vector interaction did not contribute to 
the superconducting gap equation nor to the mass gap equation unlike the case here. We have compared the pressure 
in both the cases in this regime and it appears that in this small window, existence of both the condensates has higher 
pressure than having only the quark antiquark condensates in the ground state. Beyond this window chiral symmetry 
is restored and the gap increases monotonically with chemical potential till the effect of cut off is felt and then it 
decreases. 
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FIG. 1: Masses of the quarks as a function of baryon chemical potential. Solid curve refers to masses of u and d quarks while 
the dashed curve refers to strange quark mass. 

We next discuss the results when charge neutrality conditions are imposed. As earlier, let us focus our attention 
first to the case without the diquark condensates. The results are shown in Fig. 3. Here, the d- quark mass vanishes 
earlier than that of u quarks as the baryon chemical potential is increased. The reason is that to maintain electrical 
charge neutrality conditions the d quark densities are almost twice that of u quarks which makes the d quark antiquark 
condensates to vanish. In contrast to the rather sharp fall of the masses as compared to without imposition of electrical 
neutrality conditions, the u quarks remain massive much after d quark become massless (about 80 MeV in the quark 
chemical potential window). The magnitude of electric chemical potential increases with densities till strange quarks 
begin to appear beyond which it starts decreasing so that charge neutrality conditions are maintained. 

Next we show the results with diquark condensates with neutrality conditions. For a given quark chemical potential 
/X, the interacting chemical potentials are determined using Eq. l|37|l with a trial value of (jle and /its- For high baryon 
chemical potential when chiral symmetry is restored for light quarks, the thermodynamic potential is varied with 
respect to the strange quark constituent mass after solving the superconducting gap equation Eq. l|44|l . The values 
of he and /!$ are varied so that the charge neutrality conditions Ea. (|48|l and Ea. (|49|l are satisfied. The results are 
shown in Fig. 4. The behaviour of the superconducting gap is similar to that when charge neutrality condition are not 
imposed except that the magnitude of the gap decreases. The pressure of color and electrical neutral superconducting 
phase is lower than when neutrality conditions are not imposed as shown in Fig|SJ Further, the pressure of color and 
electric neutral normal quark matter is always lower than that of the color and electrically neutral superconducting 
matter. 

At high densities when chiral symmetry is restored for the light quarks, S € — and S = —S v = — {g 2 /8)(p u — 

p d )). For S v < 0, u>i, the quasi particle energy for u quarks is always positive and hence only the d quark contribution 
will be there in the second term of the thermodynamic potential given in Ea. l|45() . Further, it is easy to show 
that such a contribution arises when the magnitude of S Ul which is half the difference of the chemical potentials of 
the quarks which condense, is the same or greater than the superconducting gap, A. When \S U \ — A, the mode 
L02 = ^/A 2 + (e 2 — v 2 ) — S, becomes gapless at the fermi sphere. For \S V \ > A, the gapless modes occur at momenta 
higher than the (average) fermi momenta. In the present calculation this is the case for baryon chemical potential 
below 1600 MeV. The occurrence of such gapless superconducting modes in neutral quark matter was first emphasised 
in Ref. |2^] . Because of this the number densities of u and d quarks participating in the superconducting phase are not 
the same in this region. This is plotted in Fig. El It may be noted that even in the gapless superconducting mode, 
the density of strange quarks is small but nonzero. 
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FIG. 2: Masses of the quarks and the superconducting gap as functions of baryon chemical potential for the case fiE =0=fie- 
Solid, dashed and dot-dashed curves refer respectively to masses of up (down) quarks, strange quark and to the superconducting 
gap. 



5. SUMMARY 



We have analysed here in a current current point interaction model, the structure of vacuum in terms of quark 
antiquark as well as diquark pairs. The methodology uses an explicit variational construct of the trial state and 
is not based on a mean field calculations. Because of the point interaction structure we could solve for the gap 
functions explicitly. The distribution functions are also determined variationally. Because of the vector interaction, 
the chemical potentials are interaction dependent and are calculated self consistently. We find that there is a small 
window in baryon chemical potential (about 80 MeV) when both the condensates are nonzero. 

To consider neutron star matter, we have imposed the condition of color and electric charge neutrality conditions 
through introduction of appropriate chemical potentials. It has been noted and emphasized earlier that projecting 
out the color singlet state from color neutral state costs negligible free energy for large enough chunk of color neutral 
matter [l3| . 

The gap reduces when color and electric neutrality conditions are imposed. The pressure with the gap is always 
higher than free quarks when charge neutrality conditions are imposed. For slightly lower densities, but large enough 
to be in the chiral symmetry restored phase, there appears to be gapless modes available when the condensing quarks 
have a difference in chemical potentials which is larger than twice the superconducting gap. In all these calculations 
we have also included the effect of self consistently determined mass for the strange quark. In fact, in the gapless 
superconducting phase, the number densities of strange quarks is also nonzero. 

We have focussed our attention here to the superconducting phase. The variational method adopted can be directly 
generalised to include color flavor locked phase and one can then make a free energy comparison regarding possibility 
of which phase would be thermodynamicaly favourable at what density. This will be particularly interesting for 
cooling of neutron stars with a CFL core. We have considered here homogeneous phase of matter. However, we can 
also consider mixed phases of matter with matter being neutral on the average. Some of these problems are being 
investigated and will be reported elsewhere [2(j . 
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FIG. 3: Masses of the quarks without superconducting phase when electrical charge neutrality condition is imposed. Solid, 
dot-dashed and dashed curves refer respectively to masses of up, down and strange quark mass. 




FIG. 4: Superconducting gap when color and electrical charge neutrality condition is imposed (solid curve). The dashed curve 
corresponds to when this condition is not imposed (i.e.^,E=0). 
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FIG. 5: Pressure as a function of baryon chemical potential when color and electrical charge neutrality condition is imposed 
(solid curve). The dashed curve corresponds to when this condition is not imposed. 
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